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The  velocity  data .obtataed tom 

by  photons.  The  velocity  vector  is  determined  from  the  displacement  of  smdl  «^ons  y 

*  6d  SZ L“^of  PW  To  taL  advantage  of  standard  data 

technique  may  be  viewed  as  aemctecula^  ^  ^  onto  a  regular  grfd  with  a  uniform 

processing  tec  ^qu ^ examine  the  ^curacy  and  noise  issues  related  to  the  use  of  various  low 

Sr  poIy™Ll  least-square  fits  for  remapping  and  the  subsequent  computation  of 
Se  datr^Tfre  infotmation  obtained  has  relevance  to  PIV  data  processing  as  well.  As  noted  by 
Sdinl  andRienot  (1993)  the  best  estimate  of  the  location  of  the  velocity  vector  acquired  through 
St of  as  PIV,  is  at  the  midpoint  of  the  displacement  vector.  Thus,  unless 

^nerifil  care  is  taken  PIV  data  are  also  initially  obtained  on  an  irregular  gnd. 

In  the  past  various  methods  have  been  used  for  remapping  randomly  spaced  velocity  data  on  0 
a  regular  grid  Among  them  are  an  inverse  distance  approach  and  a 
by  Spedding  and  Rignot  (1993).  In  this  study,  we  consider  the  use  of  2  ,3  ,or4 
to  remap  the  velocity  field  by  performing  a  local  least-squares  fit  of  fie 

within  region  of  radius  R.  This  choice  was  selected  based  on  the  work  of  Agui  and  Jimenez  (198  ), 
who  report  that  low  order  polynomial  fits  and  “kriging”  techniques  produce  the ,mos^ur^ 
representation  of  the  actual  velocity  field.  However,  no  quantitative  information  on  toe  P«rf^nce 
of  these  methods  is  presented.  Four  approaches  are  assessed  here  for  computing  the  out-of-plane 
vorticity  field  from  the  in-plane  velocity  measurements:  direct  differentiation  of  ^  P0^0™*1 
n,ed  in  the  remapping  process,  1st  and  2fld  order  finite  difference  techmques  (2  and  4  order 
accurate,  respectively)  and  an  8-point  circulation  method  on  the  regular  data.  It  should  be >  noted It  at 
the  8-point  circulation  method  is  identical  to  the  1st  order  finite  difference  method  after  an 

aPPT“:rf1T0cSned  the  uccumcy  of  the  varfous  means  to  compute  vorticity  from 
regularly  spaced  data.  Spedding  and  Rignot  (1993)  used  a  1-  order  finite  difference  technique  for 
the  inverse  distance  method  or  directly  differentiating  the  global  basis  function.  Results  indi 
that  the  Global  basis  function  produced  generally  superior  results,  however,  results  were  found  to  b 
highly  dependent  upon  the  ratio,  of  a  characteristic  length  scale,  of  the  flow,  L,  to  the  mean  spacing 
between  measurements,  8.  Abrahamson  and  Lonnes  (1995)  found  that  ccjtcu  a .on 
resulted  in  slightly  more  accurate  vorticity  results  than  using  a  least-squares  fit  to  a  model  velo  y 
field.  Luff  t*  dl(1999)  compared  the  1«  and  2“>  order  finite 

circulation  method  in  the  calculation  of  vorticity  in  the  presence  of  both  noise  and  missing  < 
points.  In  terms  of  only  the  computed  vorticity  RMS,  the  1st  order  finite  difference  technique 

produced  the  best  results.  .  , 

One  short-coming  of  the  above  mentioned  studies  is  that  only  the  random  component  ot  the 

error  field  is  examined.  Fouras  and  Soria  (1998)  found  that  the  error  in  the  vorticity  field  could  be 
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.  •  c  :t : _  divided  into  two  portions:  a  mean  bias  error  due  to  spatial  filtering,  and  a 

k  sh  mtssss = s 

S  tetSL  on  differentiating  the  fit  were  sensitive  to  the  number  of  points  used  in 
the  fit.  This  work  was  based  entirely  on  regularly  sampled  velocity  data,  issues  conn 

JSSTS3SZ,  we  include  the  effect  of  remapping  of  the  velocity  field  on  the  accuracy  of 
vorticitv  estimation  The  differentiation  of  the  polynomial  fit,  used  in  remapping,  as  a  means  of 
esdmatina^e'voScity  iTaiso  considered.  A  simulation  of  an  Oseen  vortex,  utilized  by  the  previous 
studiesTis  also^mployed  here  as  die  basis  for  comparison.  The  effect  °^uncert^^^ 
measurements  in  sLulated  by  adding  a  random  amount  of  no.se  to  each  velocity  compo  . 

Our  simulations  show  that  the  accuracy  of  the  remapping  process  depend  on  both  1/8 i  and 
norm^zM  size  of  the  region  used  in  the  least-squares  fit,  »8.  The  effect  of  mcreasmg  the  value 
of  R/S  on  the  velocity  and  vortkity  flelds  is  to  decrease  the  random  error,  >=«  “OT^e  tte  nwanb 
error.  Generally,  the  decrease  in  random  error  is  smaller  then  die  increase  m  the 
the  different  polynomial  orders  tested  produced  accurate  results  in  the  remapping  of  the  velo  ty 
for  “chotae  of  1/6  and  «i  .  For  example,  with  i/5  >  4.5  both  the  velocity  mean  bias  error 

and  random  error  are  less  than  1%  of  the  peak  velocity.  rii  fl 

Results  show  that  the  most  accurate  vorticity  results  are  achieved  by  either ^ecfly 
differentiating  the  3"  order  polynomial  fit  to  the  original  inegular  data  or  by  *e  2 .  order  f  m 
difference  technique.  Note  that  this  conclusion  is  different  fromprevious  stnd.es  because  they  either 
did  no.  consider  the  2“*  order  finite  difference  approach,  or  the  performance  was  based  only  on  the 
random  error.  Fori/8  >  4.5,  we  find  the  mean  bias  error  and  random  OTOr  m  ^  ‘me 

than  3%  and  2%,  respectively,  of  the  peak  vorticity.  A  decrease  in  i/8  causes  a  large  increase  in 
vorticity  bias  error,  error  values  higher  than  1 8%  of  the  peak  vorticity  are  found  for  i/B  -  2.5 . 
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INTRODUCTION 

Often  cracks  in  rubbery  particulate  composites 
are  found  in  the  propellant  grain  of  rocket  motors  and 
experience  mixed-mode  loading.  These  composites  are 
made  from  a  rubbery  matrix  with  a  high  volume  of  rigid 
particles  (70-80%).  The  cracks  threaten  structural 
integrity  and  can  grow  to  catastrophic  failure; 
Refinements  in  our  predictive  abilities  yield  cost 
savings  by  improving  service  life  predictions.  This 
work  uses  an  approximate  analysis  to  examine  the  crack 
behavior.  This  approach  is  practical  and  useful  because 
it  can  be  put  to  widespread  use  in  industry. 
EXPERIMENTAL  PROCEDURE 

We  used  center  cracked  specimens,  101.6 
mmx  101.6  mm  x  5.08  mm,  with  25.4  mm  center 
cracks  oriented  at  an  angle  p  with  the  horizontal. 
Different  crack  orientations  provided  specimens  with 
different  mode  mixities.  For  these  specimens,  the  phase 
angle  of  the  complex  stress  intensity  factor  was 
approximately  equal  to  the  crack  orientation  angle  /?. 

We  tested  two  specimens  at  each  of  the  crack 
orientations  (0°,  15°,  30°,  45°,  and  60°),  applying 
uniform  vertical  displacements  to  the  horizontal 
specimen  edges  at  a  constant  5 .08  mm/min  rate.  We 
measured  the  initiation  loads  and  growth  rates  from 
videotape;  we  measured  the  initial  growth  directions 
from  the  fractured  specimens  after  testing.  The 
initiation  loads  were  used  with  finite  elements  to 
determine  the  components  (KIC  and  Knc, )  of  the  initiation 
toughness. 

RESULTS 

We  obtained  all  the  results  shown  by  applying 
linear  elastic  fracture  mechanics  concepts  to  a  rubbery 
composite  that  has  viscoelastic  properties  capable  of 
high  elongations.  This  approach  is  justified  because  the 
fracture  initiation  and  growth  occur  on  the  linear  part  of 
the  stress-strain  curve  and  because  moire  experiments 
show  linear  elastic  behavior  near  cracks  in  these 
composite  specimens.  However,  because  of  the  time 
dependent  constitutive  behavior,  the  results  must  be 
applied  to  structures  with  similar  nominal  strain  rates. 
This  use  of  linear  elastic  fracture  mechanics  is  practical 
because  by  greatly  simplifying  analyses  that  would 
otherwise  be  cumbersome  and  unfeasible  it  allows  for 
widespread  use. 


The  fracture  toughness  locus  for  the 
mixed-mode  tests,  shown  in  Figure  1,  has  an  elliptical 
curve  fit  of  the  form: 
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Figure  1  Elliptical  failure  locus 

We  treated  the  parameters  KIC  and  Kac  as 
unknowns  and  determined  them  using  least-squares. 
Equation  (1)  is  regarded  as  linear  in  the  unknowns 
(l/Krc)2  and  (1  /KIIC)2,  and  the  least-squares  method  is 
applied  using  matrix  algebra  1 .  Researchers  have  used 
an  elliptical  curve  fit  previously  for  isotropic  materials. 
Note  that  the  pure  mode  II  fracture  toughness  is  smaller 
than  its  mode  I  counterpart.  Previous  explanations  that 
examined  the  fracture  toughness  locus  in  relation  to 
micromechanisms  were  based  on  observations  from 
metals,  and  do  not  seem  to  apply  to  rubbery  composites. 
Determining  the  connection  to  microstructure  is  left  for 
future  work,  but  two  possibly  relevant  phenomena  are 
suggested:  the  presence  of  rubbery  ligaments  behind  the 
crack  tip  and  the  growth  of  voids  from  particles  near  the 
crack  tip. 

We  used  the  experimentally  determined  loads 
with  finite  elements  to  determine  the  J  integral  and  then 
to  derive  KIC  and  Knc-  Using  the  crack  face 
displacement  data,  it  was  determined  that  for  these 
specimens  the  complex  stress  intensity  factor  phase 
angles  were  approximately  equal  to  the  crack 
orientation  angles.  Then,  using  Jc  =  Kc2 IE,  where  Kc  is 
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the  magnitude  of  the  complex  fracture  toughness,  and 
XF=  /?,  where  TGs  the  phase  angle,  allows  K;C  and  Knc 
to  be  calculated. 

The  initial  growth  directions,  or  kink  angles, 
were  determined  by  experiment  and  are  shown  in 
Figure  2.  This  figure  also  shows  predicted  kink  angles 
from  several  theories,  all  of  which  gave  similar 
results 2,3 .  The  large  deformations  during  actual  loading 
made  it  difficult  to  unambiguously  determine  the  kink 
angles,  so  we  measured  them  after  testing  by  projecting 
the  fractured  surfaces  onto  a  large  screen. 


a  max  tensile  stress 
«  max  energy  release  rate 
i  max  strain  energy  density 
*  propellant  experiments 


Figure  2  Kink  angles  vs.  crack  orientation  angle 


We  measured  the  growth  rates  from  videotape 
using  the  approach  of  Figure  3 .  During  most  of  the 
growth,  the  angled  crack  of  Figure  3a  experiences 
mode  I  growth  that  can  be  characterized  by  the  effective 
crack  length  and  related  mode  I  stress  intensity  factor  of 
Figure  3b.  The  crack  growth  rate  and  Ki  can  be  related 
through  a  power  law  relation  such  as: 
da  „ 

— 4-  =  CK 7  (2) 


<») 


This  gives  nearly  identical  results  for  all  the 
tested  mode  mixities.  Because  of  this,  the  growth  of  a 
mixed-mode  crack  with  a  geometry  shown  in  Figure  4a 
can  be,  modeled  using  the  simplified  geometry  of 
Figure  4b.  Figure  5  shows  the  aggregate  data  and 
resulting  curve  fit;  the  parameters  in  eqn  (2)  are 
C  =  1.79  x  10‘6  and  m  =  2.73.  Results  for  a  mode  I 


specimen  gave  the  nearly  identical  results  of 
C  =  1.85  x  1CT6  and  m  =  2.74. 


Figure  5  Effective  crack  growth  rate  vs.  Ki 

CONCLUSIONS 

Although  rubbery  particulate  composites  have 
viscoelastic  properties,  high  elongations,  and 
complicated  failure  mechanisms,  they  can  be  studied, 
for  a  given  nominal  strain  rate,  using  the  principles  of 
linear  elastic  fracture  mechanics.  When  analyzed  like 
this,  the  fracture  locus  is  elliptical.  The  initial  crack 
growth  angles  match  the  strain  energy  density 
predictions  best,  although  for  all  but  the  highest  mode 
mixities,  other  theories  gave  nearly  identical  results. 

The  crack  growth  rates  can  be  predicted  using  an 
approximate  mode  I  approach.  The  result  is  a  simple 
approach  that  circumvents  more  detailed  and  difficult 
analyses.  Future  work  could  include  experimental 
methods  to  study  the  causes  of  fracture  in  these 
composites,  especially  the  micromechanisms  near  the 
crack  tip  at  various  mode  mixities. 
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Figure  3  Modeling  mixed  mode  crack  growth 


